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Abstract 

The ageing algebra is a local dynamical symmetry of many ageing systems, far from 
equilibrium, and with a dynamical exponent z = 2. Here, new representations for an 
integer dynamical exponent z = n are constructed, which act non-locally on the physical 
scaling operators. The new mathematical mechanism which makes the infinitesimal gener- 
ators of the ageing algebra dynamical symmetries, is explicitly discussed for a n-dependent 
family of linear equations of motion for the order-parameter. Finite transformations are 
derived through the exponentiation of the infinitesimal generators and it is proposed to 
interpret them in terms of the transformation of distributions of spatio-temporal coor- 
dinates. The two-point functions which transform co-variantly under the new represen- 
tations are computed, which quite distinct forms for n even and n odd. Depending on 
the sign of the dimensionful mass parameter, the two-point scaling functions either decay 
monotonously or in an oscillatory way towards zero. 



1 Introduction 



Non-relativistic space-time transformations have recently met with a lot of interest - in addi- 
tion to fields such as hydrodynamics [4"3~1 152] , they have been playing an increasing role in the 
analysis of the long-time behaviour of strongly interacting many-body systems far from equi- 
librium [12} |27] and even more recently have arisen in non-relativistic limits of the AdS / CFT 
correspondence, see e.g. [U EH El [32] , with interesting applications to cold atoms [HI [15]. One 
of the ingredients for physically interesting sets of space-time transformations appears to be 
some kind of conformal invariance and recently, a classification of non-relativistic conformal 
space-time transformations was presented [H]. An important sub-set of these is characterised 
by a finite dynamical exponent < z < oo. Indeed, the list of sets of admissible generators 
which close into a Lie algebra is a rather short one, in d + 1 space-time dimensions: 

1. the conformal algebra itself, in al + 1 dimensions, with z = 1. 

2. the Schrddinger algebra, which was discovered by Lie in 1881 as a dynamical symmetry 
algebra of the ID free diffusion equation^ The dynamical exponent is z = 2. This algebra 
was rediscovered in physics as a symmetry of free non-relativistic particles several times 
around 1970, including [331 [T8| HUJ EI]- Non-linear examples of Schrddinger-invariant 
equations include the Navier-Stokes equation [131 EH HZ] or Burger's equation [HI 130] . 
see e.g. [27] and refs. therein for a short summary of further examples. Schrddinger 
symmetry (or rather its sub-algebra with time-translations left out) also arises in the 
far-from-equilibrium dynamics of statistical systems [21] , for example in simple magnets 
quenched to a temperature T < T c below the critical temperature T c > from a fully 
disordered initial state when z = 2 is known [HI [T21 EZ] • 

3. the conformal galilean algebra CGA(g?) appears to have been first identified in [20J, but 
was independently rediscovered in different contexts [22, 39J. It is usually obtained, by 
a contraction, as the non-relativistic limit of the (d + 2)-dimensional conformal algebra 
(itself obtained by a non-relativistic holographic construction) [2U E7J [H [21 EH ESI E2] • 
In d = 1 space dimension, there exists an infinite-dimensional extension which can be 
constructed from a contraction of a pair of commuting Virasoro algebras [231 EHl 12] . In 
most representations, one has z = 1, but representations with z = 2 are also known [26] . 

4. in d = 2 space dimensions, there exists the exotic conformal galilean algebra ECGA, which 
is the central extension of the non-semi-simple CGA(2) [36]. Although one may readily 
identify linear equations invariant under ECGA [37], the construction of invariant non- 
linear equations is not straightforward [521 ELI]- The dynamical exponent is z — 1. See 
[2H] for a recent review and the relationship with non-commutative mechanics. 

5. finally, there exists for d = 1 a closed algebra with z = | [23]- It is not yet clear how 
this might fit into the general scheme of [TS] , since it does not contain the full conformal 
structure and furthermore its generators contain fractional space derivatives. 

This short list illustrates the practical difficulty of constructing sets of 'conformal' space-time 
transformations for a generic dynamical exponent z ^ 1, 2 such that a Lie algebra is obtained. 

1 Jacobi already wrote down en passant in 1842/43 the elements of the corresponding Lie group as symmetries 
of the Hamilton- Jacobi equation of a particle with an inverse-square potential, see |20j . 
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It is at present not fully understood how to construct a dynamical symmetry even for a simple 
linear equation of the form (where 1,2) 



Stl>{t, r) := (zfidt - d*)i/>(t, r) = (1.1) 

which arises as one of the most simple equations of motion of the order-parameter in studies 
of ageing far from equilibrium (Till E]- Indeed, current attempts to find further dynamical 
symmetries of eq. (11 .11) beyond the obvious translation-, dilatation- and rotation-symmetries 
(if d > 1) only succeed at the price that the further generators must be required to vanish on 
certain states (which are then declared to be the 'physical' ones) [22l EHJ, [27] . Furthermore, these 
generators cannot, in general, be expressed in terms of first-order differential operators ('vector 
fields' in mathematical terminology). In the context of statistical physics, the order- parameter 
does not really satisfy a deterministic equation, but rather the r.h.s of (11.11) is replaced by a 
random noise term, which leads to a Langevin equation. However, since the non-relativistic 
algebras mentioned above are all non-semi-simple and their representations are projective, it 
is possible to study first the symmetries of the deterministic equation (11. ip and then use the 
resulting Bargman super-selection rules [5] in order to reduce the calculation of any average to 
the calculation of averages within the deterministic part of the theory as defined by (II. ip [44J. 
This procedure works not only for thermal noises and a simple diffusion equation with z = 2, 
but can be generalised to generic values of z and fairly general noises, such as they may arise 
in reaction-diffusion systems |6j EHJ [7J El [13] , see [27] for a systematic presentation. 

In this paper, we shall explore properties of a new kind of representations of the common sub- 
algebra agz(d) of the Schrodinger and conformal galilean algebras. For the sake of notational 
simplicity, we shall restrict from now on to d = 1 spatial dimensions. Then the standard 
representation, on sufficiently different iable space-time functions f(t,r), of the Lie algebra 
agc(l) := (X 0A ,Y ± x, M ) is given by 

I x 

X = —td t — -rd r — — dilatation 

X\ = — 1 2 <9( — tr<9 r — — r 2 — (a; + £)t special transformation 

Y_x = —d r space-translations (1.2) 

Yx = —td r — Air Galilei-transformation 

2 

M = —M phase shift 



See (12. 2p for the commutators. This representation is characterised by the 'mass' M. and the 
pair of scaling dimensions (x,£) whose values depend on the scaling operator on which these 
generators act. If one defines the Schrodinger operator 

S := 2Md t -d 2 r +2M (x + £- j (1.3) 

then the equation Si/)(t, r) = has oge(l) as dynamical symmetry, because of the commutation 
relations 

[S,Y±x] = [S, M ] = , [S,X ] = -5 , [S,X 1 ] = -2tS (1.4) 

which imply that any solution of Sip = is mapped onto another solution of the same equation. 
A physical example for ( 11 .3p is given by the relaxation kinetics of the spherical model, or equiva- 
lently the N — > oo limit of the O(N) model, after a quench to a temperature T < T c at or below 
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its critical temperature T c > [TH]. The representation ( II. 2p has a dynamical exponent z = 2 
and acts locally on the space-time coordinates. We point out that because time-translations 
(with a generator X_i = —d t ) are not included and hence a system with an age-symmetry is 
not at a stationary state, the scaling dimension £ arises as a further universal characteristics 
of the relaxation process. By the transformation <f)(t,r) = t^$(i, r), the physical age-quasi- 
primary field with scaling dimensions (x^, can be related to the Schrodinger-quasi-primary 
field $ with scaling dimensions (x$ = x$ + 2^, = 0) [25] and in the transformed eq. (II. 3p 
x + £ is replaced by If this option is chosen, one must give up the identity of physical and 
quasi-primary scaling operators, familiar from conformal invariance, which holds for station- 
ary, equilibrium systems. The potential term in ( 11. 3ft can of course be eliminated by a similar 
transformation or else by imposing the constraint x + £ = |. 

When trying to extend ( 11.2ft to a representation of cga(1) D age(l), the extra generators 
are not necessarily first-order differential operators [24J. We believe that this fact should be 
taken seriously and its consequences studied. For this reason, and in order to find further 
representations of age(l) with different values of z, we shall give in section 2 non-local repre- 
sentations of age(l), which admit any integer value z = n e N, but which cannot be reduced 
to first-order differential operators. As we shall see, closure of these representations requires 
to restrict the representation space by performing a quotient with respect to the Schrodinger 
equation Sip = 0. In section 3, we address the question how to interpret geometrically such in- 
finitesimal generators by explicitly constructing the finite space-time transformations given by 
the exponentiated generators. The examples studied here suggest that these finite transforma- 
tions might be viewed as transformations of distributions of space-time coordinates, instead of 
precise transformations of the coordinates. An appendix compares this with the finite Galilei- 
transformations. Next, in section 4, we derive the co- variant two-point functions which depend 
strongly on the parity of n. We conclude in section 5. 



2 Non-local representation of the ageing algebra age(l) 

Consider a dynamical exponent with integer values 2 < z = n e N. The generators of age(l) 
we are interested in take the form 

An = to t ret 

u 2 2 2 

7i i 

X 1 = --t 2 d t d?~ 2 -trd?- 1 - -fir 2 -(x + Otd?- 2 

Y_i = -d r (2.1) 

Yi = -td^-fir 
Mo = -ii 

and satisfy the commutators of age(l), of which we give the non- vanishing ones 

[A ,F ± i] = t\y ± . , [X Q ,X 1 ] = -X, , [Yk,y_i] = M . (2.2) 
However, there is a notable exception, namely 

n — 2 

[X u Y k ] = —t 2 dr 3 S (2.3) 
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with the Schrodinger operator 



The generators (12. ip form a dynamical symmetry of the Schrodinger equation Sip(t,r) = 0, as 
can be seen from the commutators 

[S,Y ± i] = [S, M ] = , [S,X ] = --S , [S,X 1 ] = -ntVr 2 S. (2.5) 
2 z 

In order to close the representation (12.11) . we must restrict the function space modulo solutions 
of Sip = 00 Then a natural function space for our purposes is F := C 1 (R + , C n (R))/ ~, the 
space of functions which are continuously differentiable in time and n times differentiable in 
space or alternatively F = C 1 (R + , H n (R))/ ~, where the Sobolev space H n (R) of n-times 
differentiable functions which together with their derivatives are also square-integrable is used 
such that Fourier transforms with respect to r exist; and at the end with the quotient taken 
with respect to the Schrodinger equation Sip = 0. 

Restricted to the space F, the generators (12. ip give for each integer n > 2 a non-local 
representation of Ofle(l) which is a dynamical symmetry of the Schrodinger equation Sip = 0. 

3 Finite transformations 

Besides the usual local generators of dilatations X , of spatial translations Y_i and of phase 
shifts Mo, the representation (12. ip also contains the non-local generators Yi,X\ whose effect 
cannot be interpreted as a simple space-time coordinate transformation t h- >■ t'(t, r), r i— > r'(t, r). 
On the other hand, we can still write the formal Lie series F(e,t,r) = e~ eYl / 2 F(0,t,r) and 
F(e,t,r) = e~ eXl F(0,t,r). They are given as the solutions of the two initial-value problems 

{de-td 1 ^- 1 -fiAF(e,t,r) = , F(0, t, r) = 0(t, r) (3.1) 
(d E - ^t 2 d t d?- 2 ~ trd?- 1 - xtd?~ 2 - ^/xr 2 ) F(e, t,r) = , F(0, t, r) = <f>(t, r) (3.2) 

such that the initial function <p £ F . 

In tables [T] and (2J we illustrate these Lie series for the choices 4>(t,r) = t m and 4>{t,r) = r k 
with m £ N and 1 < k < n — 1, which for /i = solve the Schrodinger equation <S0(i, r) = 0. 
Comparison with the effects of the local Galilei- and special Schrodinger transformation shows 
important differences. For example, although the spatial coordinate r is left invariant by both 
generators when n > 2, this does not imply that these generators would not generate any 
spatial transformation, as we see from the transformation behaviour of the higher powers of r. 
While in the local case n = 2, the transformation of the powers r k is simply given by taking 
the corresponding power of the transformation law of r itself, this is no longer true in the 
non-local cases n > 2. While the action of the generators Yi and X\, in our example, cannot 
be interpreted in terms of a local coordinate transformation, the results look reminiscent to a 

2 This can be done by the following construction: to functions / and g are said to be equivalent, written 
/ ~ g, if there is a sufficiently differentiable function A(t, r) such that f(t, r) = g(t, r) +A(t, r) and SA(t, r) = 0. 
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Table 1: Comparison of the finite transformations e~ eYl / 2 (fi(t, r) for the generalised, non-local 
Galilei-transformation when z = n > 2 with the standard local Galilei-transformation for 
z = n = 2. The initial distribution G J 7 and \x = 0. 



<j>{t,r) 


non-local, n > 2 local, n = 2 






r k (r + te) k 
r n-i + ( n _l)!te (r + te)™" 1 


m G N 

1 < jfe < 71-2 



Table 2: Comparison of the finite transformations e~ eXl <j)(t,r) for the generalised, non-local 
special Schrodinger-transformation when z = n = 3 or 4 with the standard local special 
Schrodinger-transformation for z = n = 2. The initial distribution G J 7 and /i = 0. 



4>(t,r) 


non-local 

n = 3 n = 4 


local 

n = 2 




r 
r 2 

^3 


r + [x + ^)te r 
r 2 + 2(x + ^ + l)tre r 2 + 2(x + £)te 
+i(x + e + l)(2x + 2e + 3)t 2 e 2 

r 3 + 6(x + £ + l)tre 


t m /(l -te) m+x+ Z 

r/(l-te) 1+x+ t 

r/(l-te) 2+x+ t 

r/(l - te) 3+a;+ « 


m G N 



transformation of a statistical distribution, where the first moment happens to be invariant, but 
the higher ones change. Therefore, these examples suggest that a better interpretation might 
be to consider a transformation of an initial distribution of spatial (or temporal) coordinates, 
where (f>(t,r) would then take the role of a distribution function. 

In what follows, we give further results on the transformation of <ft(t, r) and discuss possible 
consequences for an interpretation. In order to keep the expressions to a manageable size, 
we shall concentrate on the two cases z = 3 and z — 4. These are the values of z in the 
Bray-Rutenberg theory of the growth of the relevant time-dependent length scale L(t) ~ t l l z 
in 0(n)-symmetric systems with a conserved order parameter and quenched to T <T C [9]. 

3.1 The case z = n = 3 

We now give the full transformation laws of the distribution <p(t,r). We begin with the gener- 
alised Galilei transformation (13.11) and use the Fourier representation 

F(e, t, r) = —L= [ dk e ikr F(e, t, k) (3.3) 

V 27T Jm. 

This leads to the equation (d £ +tk 2 —ifidk)F = 0. Letting v := e— ik/fi and F(e, t, k) = G(v,t,k), 
it readily follows that G(v,t,k) = G (v , t) exp(— it/c 3 /(3/i)) and where G (v,t) must be found 
from the initial condition, with the result Go(v,t) = <p(t, ifik) exp(t/i 2 /c 3 /3). This gives the 
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transformed distribution in Fourier space 



4>(t, k) i — y F(e, t, k) = <f>(t, k + ifie) exp 



-tk 2 e - ifitke 2 + -fi 2 te 3 
3 



(3.4) 



and finally in direct space, after having performed the integral over k, the general solution to 
eq. (13. ip becomes 



F(e, t, r) = 1 / dr' (f>(t, r') exp — — ( (r - r' - tfie 2 ) 2 - Afxtr'e 2 - -/i 2 t V ] 
\fA-nte Jr I 4te \ 3 J 



(3.5) 



Setting ji = 0, we obtain the entries in table [TJ Up to the /x-dependent terms, eq. (13 .5p is a 
convolution of the initial distribution with a gaussian and using the form (13.41) . it is readily 
checked that the group property holds true. 

Specifically, we list some examples of finite transformations when [i ^ 0: 

r' = (r + /ite 2 )e^ r+Me3/3 
(r 2 )' = ((r + /ie 2 t) 2 + 2te)e^ + ^ 3 / 3 

(t + ^r 3 J = ^+|(r + /ie 2 t) 3 + 3/ite(r + /ite 2 ))e^ r+ ^ 3/3 . (3.6) 

We have checked that these solve (13. ip . as well as the Schrodinger equation with n = 3. 

In particular, if one tentatively interprets 0(r) as a probability distribution such that 
J R dr(f)(t, r) = 1, this normalisation condition remains unchanged for /i = 0, viz. J R dr F(e, t, r) |^_ Q 
1. Furthermore, one may consider 

0(t, k) = (e- ikr ) = —L= [ dr e" ifcr 0(t, r) (3.7) 



/2tt 

as the associated characteristic function. For example, if we consider a shifted gaussian with 
characteristic function <f)(t, k) = exp(— Xk 2 + ijk), this transforms into 

?(t, k) ^ F(e, t, k) = e -( A + te ) fc2+i ^ e -ik(^Xe+^) e ^(X + te/3)e 2 -^e ^ 

For fi = the centre stays unchanged at 7, while the width becomes A y A + te. Gaussian 
distributions are therefore co-variant under the generalised Galilei generator Yi with fi = 0. 
However, since the gaussian distribution is not a solution of the Schrodinger equation with 
ii / 2, one can realise a gaussian distribution at best as an initial condition which has to be 
evolved in time. This illustrates the non-trivial constraint of remaining within the reduced 
function space, introduced in section 2. 

The integration of the generalised special transformation (13. 2 p runs along similar lines. For 
the sake of brevity, we set /1 = from now on. In Fourier space, we introduce the new variables 
u:=e + 2\/{kt) and v := kt 2 ' 3 and find F{e, t, k) = G{t, u, v) = t 2 ^- x -^/ 3 G (u, v). The as yet 
undetermined function Go is related to the initial distribution via 

G (a, P) = -±= (J^j 2(XH ^ fdr' exp (y'a 2 /) ( ( — Q ) , r' ) (3.9) 
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Hence, the final form for the solution of (13.21) with // = reads 



F(e, 



t, r ) = -L I dkdr'e ik ^ (l + ^ 



2tt 



2i 



2(1-2-0 



tke\ 



In particular, the entries in table [2] are recovered. 

When we consider a gaussian distribution, we find the formal transformation 



(3.10) 



</>(t,k) 



*->F(e,t,k) = 1 + 



~2T 



2(l-x-0 



7T 



-\ eB (tk)k 2 



(3.11; 



but now with a fc-dependent effective width A e ff(tfc) = A(l + tfe/(2i)) 4 . In contrast to the 
generalised Galilei transformation studied before, the transformation law also depends on the 
value of the scaling dimension x + £ and we see that the companion factor reduces to unity, 
if x + £ = 1, that is precisely when the time-dependent potential term in the Schrodinger 
equation Sip = vanishes. Again, a gaussian distribution can at best be realised as an initial 
distribution. 

Alternatively, one may implement the constraint of resting in the reduced function space 
of solutions of the Schrodinger equation directly, which we now illustrate for /i / 0. Together 
with eq. (I3.2p . we must require the Schrodinger equation 



SF(e,t,r) = ^d t + (x + ^-l)^-dA F(e,t,r) = 
This system of equations is best solved in Fourier space, where we have 



0. 



(3.12) 



3i. 



d £ - ^kt 2 d t + itk 2 d k - i{x - 2)tk + |^ ) F{e, t, k) 







(3.13) 



(3fid t + ik 3 + 2n(x + £-l)r 1 ) F(e,t,k) = 
together with the initial condition -F(0,t, k) = 0(t, k). The second of these is solved by 

F(e,t,k) = r 2 ^" 1 )/ 3 e - ifc3 */^ /(e, k) . (3. 14) 

and the first condition (13 . 13[) then leads to a diffusion equation 



d, + ^dt)f(e,k)=0. 



(3.15) 



where the diffusion constant is given by — /i/2. Standard methods give the general solution / 
and using (I3.14p we have formally 



F(e,t,k) 



F(0,t,f)exp 



2jj,e 3/i 



Going back to direct space, we finally have (using analytic continuation where necessary) 



(3.16) 



F(e,t,r) 



1 



-8ir 3 fie 



dr'dk d£(f)(t,r')e 



l\ \(k-l)r+\t{r-r') 



exp 



(k-£) 



it 



i-(k - t){k 2 + kt + f) 

2jJLE 6jJ 



1 : I dr' dm (f>(t, r' + r- tm 2 / ^yMr-t^/^+mV^s) f di exp [_ i£r > _ uHm/fx] 



y/—&7T J fl6 

1 



dm 



dr' — -— (fi(t, r' + r — tm 2 1 /z) exp 
m L ' 2 



., v \\xr' 2 m? 
im I r m H h 



t 

3/x 



Atm 2 fie 



(3.17) 



where in the last step the Fresnel integrals were used. We remark that the scaling dimension 
x + £ does not appear explicitly. 

For illustration, we write down the transformed time, which can be derived as follows. From 
(I3.17p . we have with <f)(t,r) = t, carrying out first the integral over r' via a Fresnel integral 



t' 



t 



y/-8ir 2 iet 
t 



dr' e i/f// ( 4im ) r ' ; 



dm exp 



2tt //i\V3 



4 fie V t 



exp 



it 

3/i 

2te 



dm 

m l/2 

3 1 
2ite 
1 



m 3 - 



^m(r-(t/3n)m 2 )+m 2 /(2/ze) 



2 3wr 
m m 



Ai 



(3.18) 



where Ai is the Airy function and in the second line we performed a shift in the integration 
variable in order to eliminate the terms ~ m 2 in the exponential. 



3.2 The case z = n = 4 

The finite form of the generalised Galilei transformation is found by solving (13.1 ft along the 
same lines as for the case z = 3. In Fourier space, we obtain 

$(t, k) i-> F(e, t, k) = e ^ 3e4/4 exp \it/i 2 £ s k + ~t/ie 2 k 2 - itek 3 ^j J(t, k + i/ie) (3.19) 

quite analogous to (13 .4p . From this, we find in direct space 

F(£,t,r) = eX P(^ 3£ V4) f ^ ^ [ d]( jkir-r'+t^+Zt^k'-itek^ (3 2Q) 

2tT Jr Jr 

Setting /i = 0, the results in table [T] can be recovered and we also have the same conservation 
of the normalisation, when fi = 0. Some explicit examples for transformations with fi ^ read 

(r 2 )' = ((r + /i 2 e 3 t) 2 + 3/ite 2 )e^ 3fe4 / 4 

(r 3 )' = ((r + /i 2 e 3 t) 3 + 9//te 2 (r + /i 2 te 3 ) + 6te) e ^ 3te4 / 4 . (3.21) 

Next, we integrate the special generator Xi by solving (13. 2p . Again, this is best done in 
Fourier space and we set /i = for brevity. Since with respect to the case n = 3, some 
subtleties arise, we proceed step by step. First, we introduce the new variable v = kt 1 ! 2 and 
set F(e,t, k) = G(e,t,v) which satisfies the equation 

(d £ + 2v 2 td t + (x + £- 3)v 2 )G(e, t,v) = (3.22) 

This in turn is solved by setting u = e — hit/(2y 2 ) and we find 

G(e,t,v)=g (u,v)t-^-^ 2 (3.23) 



S 



where ga{u, v) is determined from the initial condition 



g (u,v) = e~ (x+ ^ uv2 $(e- 2uv \e uv2 v 



Using the inverse transformations t = e , k = e uv v, the final form is 

~(xH~3)etk 2 2f e -2stk 2 t e etk 2 A _ 



(3.24) 



(3.25) 



<f>(t, k) \-> F(e,t, k) = e 
and this gives in direct space (with /i = 0) 

F(e,t,r) = — [ dkdr' e^-^-^' 2 ^ 2 (j)(e- 2£tk2 t, 

from which the corresponding entries in table El follow! The main difference with respect to 
( 13.1 Op is the exponential rescaling of time and space. The case /i ^ can be treated in the 
same way as n = 3 case. We omit the calculation. 



-2rffc 2 + -etk 2 I 



(3.26) 



4 Covariant two-point functions 

We now derive the form of the co- variant two-point function 

F = F{t x , t 2 ; n, r 2 ) = (^(h, n)0 2 (* a , r 2 )> (4.1) 

and where the scaling operators <pi have scaling dimension (xi, and mass /ij. The co- variance 
of F is expressed by the conditions X^F = 0, where is the two-body extension of the 
generators X e ofle(l) constructed in section 2. 

Spatial translation-invariance Y"_iF = leads to F — F(ti,t2,r), with r = r\ — r 2 . The 
mass-invariance M F = gives the Bargman super- select ion rule /ii + /i 2 = 0. The requirement 
of generalised Galilei- invariance YiF = leads to, using again the Bargman super-selection 
rule 







n-l 
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0r$ 



T " ^) ^ = (" (*i + (-l) n ^* 2 ) - Mrr) F = 

(4.2) 



It follows that one must distinguish between the cases (i) n even and (ii) n odd. 



1. n even . We rewrite the two-point function as 

F = F{u,v,r) , u:=t 1 -t 2 , v:=t 1 /t 2 , r := r 1 - r 2 (4.3) 
and obtain from the three co- variance conditions YiF — 0, XqF = and X\F = the equations 



\-iid?- 1 - mr] F = 
F = 



ra _ 1 _ xi + x 2 

uo u ro r 

2 2 2 



2 U v-1 v r 



n 



uvd v d. 



n-2 



UV 



rd n ~ l 

v — 1 



UV 

v-1 



■a 



n-2 



{x2 + &)-^—dr 2 -\lnr 2 
v — 1 2 



(4.4) 
(4.5) 



(4.6) 



! A11 entries in tables Q] and [2] can be checked by direct substitution. 
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Acting with 2 on (I4.5p . eq. ( 14. 6 p can be simplified to 

d n r ~ 2 (nvd v + — V — (xi -x 2 + 2£ 1 -n + 2) + - (x 2 - x x + 2£ 2 - n + 2) j F = (4.7) 

\ f — 1 f — 1 / 

It clear that each of the equations (I4.4ll4.5ll4.7p will fix the dependence of F — F(u,v,r) on 
one of its variables. In fact, the scaling form obtained from eqs. ( I4.4ll4.5p implies that the 
dependence of F on v factorises such that the (n — 2)-fold derivative in (14. 7p can be dropped. It 
can also be explicitly checked that the closure condition of our representation is automatically 
satisfied, as it should be. We find the following scaling form 

F(u, V, r) = t ~ ixi+X2)/n (v - [(*l+* a )/2-Ki+6-"+2] u -i[* a -* 1 +%-n+2] f ( 4 . 8 ) 

where the form of the last scaling function / = f(y) follows form eq. (14. 4p 

^r Z + ^/(^) = o ( 4 -9) 

2. n odd . We rewrite the two-point function as 

F = F(u,v,r) , u:=t l +t 2 , v\=t x jt 2 , r:=r x -r 2 (4.10) 
and now obtain from the three co-variance conditions YiF = 0, XqF = and X\F = again 

2 

the equations (I4.4|4.5p . of course with the modified relationship between u and ti j2 , while (14. 6 p 
is replaced by 



2 v + 1 2 f+1 

-(^ + so^ar 2 - + ^)^i 5 "" 2 - 

Using again ( I4.5p . we find the more simple condition 



(4.11) 



d n - 2 (nvd v + — H— (xi - x 2 + 26 - n + 2) + — ■ (ar 2 - xi + 2£ 2 - n + 2) ] F = (4.12) 
\ v + 1 v + 1 J 

This leads to the scaling form 

F(U, V, r) = t - (Xl+X2)/n (V + l)-|[(-i+-2)/2+6+6-n+2] v -l[ X2 - Xl+il -^ f ( rM -l/n) ^3) 

and where the scaling function f(y) is again given by eq. (14.91) . 

It remains to discuss the remaining scaling function f(y). The general solution of ( 14. 9 p is 

ti ^ i v 2 + £ 3 + £ n + £ /i iy n \ ... 

/ y = 22 far iFn-i i; , > • • • , ; 4.14 

where \F n _\ are generalised hyper-geometric functions and the are normalisation constants. 
On this, physically reasonable boundary conditions must be imposed, especially lim^oo f(y) = 
0. It may be more instructive, however, to look at explicit examples. 
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r/(tj+t 2 ) 



3 

1/3 



Figure 1: Scaling function f(y) in the case z = n = 3, normalised to /(0) = 1. The solid line 
gives the behaviour if fix = 1 > 0, while the broken lines indicate the behaviour, for Hi = —8 < 
and several values of k, of the function f(y) = (Ai(|/ii| 1/3 ?/) + fcBi(|/xi| 1/3 ?/))/(Ai(0)) + fcBi(O)). 



1. n = 3 . In this case, eq. ( 14. 9 p reduces essentially to Airy's equation and the solutions can be 
compactly expressed in terms of Airy's functions and the normalisation constants f± }2 



f(y) = fite{-t4 /a y) ; /ii>0 

f(y) = f 1 Ai(\ f i 1 \ 1 ^y)+f 2 Bi(\ f i 1 \ 1 / 3 y) ; ^ < 



(4.15) 



For Hi > 0, the second independent solution of (14. 9 p was suppressed, since it diverges for 
y — > oo. Figure [1] illustrates the behaviour of the scaling function, for positive and negative 
values of H\. We discuss the shape of the scaling functions below. 

2. n = 4 . The solution of ( 14. 9 p now takes the more simple form 



f(y)=f 00 F 2 ( — 



3 5 {i x y A > > 
+ Ay 0-P2 t> t; — ^r- + A 2/ 0^2 



4'4' 



64 



5 3 /iiy 4 



4'2 : 



64 



(4.16) 



This may be analysed using the leading asymptotic behaviour of the hyper-geometric function 
0-F2, which may be read off from Wright's formulae [50] 



27TV 3 



F 2 (a,b; -z) 



'3^ 



ttV3 I 2 



7T 



+ -(l-a-6)| (4.17) 



For both /ii > and /xi < 0, this implies that the function f(y) diverges exponentially fast as 
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Figure 2: Scaling function f(y) in the case z = n = 4, normalised to /(0) = 1. The thick solid 
line gives the behaviour if fii = 1 > 0. The broken lines and the grey line indicate the behaviour, 
for \i\ = — 8 < and several values of k, of the function f(y) = (J-x(y) + kjF^iy)) / '(J-i(0) + 
^(0)), with the JFi{y) defined in flTOD . 



y — > 00. We absorb this divergence by choosing the constants /q 1 2 accordingly and then find 

f(y) = fo 



0-^2 



1 3 

2' 4' 



64 



2T(3/4) 1/4 



r(i/2) 

r(3/4) 
r(i/4) 



1/2 2 77 



/3 5 




! U'4' 


64 


/5 3 




U'2' 


64 



f{y) = fo)[Fi{y) + kT 2 {y)\ 
where fo and are normalisation constants, is a free parameter and 



Hi > (4.1* 
/Ui < 



^i(v) 



/ii| 1/4 2/ 0^2 



3 5 
4'4 ; 



m\y 

64 



1/^1 



V 2 ,,, 2 



y oF2 f I'2 ; ~64~ 



r(i/2 
r(3/4) 

£(1/4 

r(3/4) 



0-^2 





3 


l/^i|2/ 4 \ 




4 ; 


64 J 



0-^2 



A 3. 




V 2' 4' 


64 J 



(4.19) 



The behaviour of these scaling functions is illustrated in figure [2J We observe that once hav- 
ing eliminated the asymptotically leading term eq. (I4.17p . the physically required boundary 
condition lim^oo f(y) = is satisfied. 

Comparing figures [1] and (2J we notice that although the scaling function satisfies for both n 
odd and n even the same differential equation (j4.9j) . the interpretation of the scaling variable 
is different. Indeed, for n = 3, the time difference t\—tz enters into y, whereas for n = 4 
it is the sum ti + 1 2 . Furthermore, we see that for fii > 0, only a single independent solution 
remains, which decreases from /(0) = 1 monotonously and very rapidly towards zero when y is 
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increased. On the other hand, for fii < 0, we find two independent admissible solutions whose 
decay towards zero is an oscillatory function of y (for n = 3, the decay should be algebraic, 
whereas it looks to be (stretched) exponential for n — 4). This feature may allow to distinguish 
at least qualitatively between two physically distinct situations with z > 2: 

• non-equilibrium relaxation kinetics with a conserved order-parameter (model B dynam- 
ics). Below the critical point, viz. T < T c , in systems with a global 0(n)-symmetry it 
is known that z = 3 for a scalar order-parameter (n = 1), and z = 4 for vector order- 
parameters (n > 2) [9J. At criticality z = A — r] = 4 — | t^W £ 2 + 0(e 3 ) in d = 4 — e 
dimensions [ST]. In these cases, scaling functions are generically seen to be oscillating^ 

• in critical dynamics, viz. T = T c , and without any conservation law on the order- 
parameter (model A dynamics), the dynamical exponent z > 2 [ST]. Here, the decay 
of the scaling functions is in general monotonous. 

Our results suggest that these physically distinct cases, even with the same value of z, might 
be distinguished through the sign of the dimensionful parameter such that /ii > repro- 
duces the monotonous decay seen in critical dynamics (model A) whereas [i\ < leads to the 
oscillatory decay found in conserved systems (model B). 



5 Conclusions 

This work has been motivated by the persistent difficulties to construct non-trivial Lie al- 
gebras of space-time transformations. We believe that the possibility of finding Lie algebra 
generators which cannot be expressed as vector fields merits serious consideration. We have 
constructed new representations of the ageing algebra age(l), corresponding to an integer dy- 
namical exponent z = n > 2 to explore the mathematical structure of dynamical symmetries 
whose infinitesimal generators are no longer described by the usual vector fields involving only 
first-order differential operators. Provided that we restrict the admissible function space to 
the solution space of the Schrodinger equation Sip = 0, and thereby somewhat relax the re- 
quirements of a dynamical symmetry, we have given an explicit n-dependent family of linear 
partial differential equations which are indeed age(l)-invariant in the sense introduced here. 
An important open question is how to extend this to non-linear equations. 

The non-local infinitesimal generators of Ofle(l) contain higher-order differential operators. 
Their exponentiation does not lead to local spatio-temporal coordinate transformations and we 
have considered the possibility that a better interpretation might be formulated in terms of 
transformation rules for distributions of spatio-temporal coordinates. Several examples of such 
transformation rules have been derived. 

Finally, we also studied the scaling form of co-variant two-point functions. Surprisingly, 
for z = n even the scaling forms are compatible with the expectations of a two-time response 
function (as it is usually the case in present theories of local scale-invariance in ageing systems) 
since they depend on the time difference t\ —t<i- On the other hand, this is not so for z = n 
odd, where the arguments of the scaling functions are much more reminiscent of co-variant 

4 For example, the scaling function ^(y) m (|4.19l) reproduces the exactly known two-time response in the 
3-D Mullins- Herring model of surface growth with a conserved order-parameter [45] . 
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two-time correlators, since they contain the sum t± + ti. We have also seen that the shape 
of the space-dependent part of the scaling functions can at least qualitatively account for the 
different forms found for non-conserved (model A) dynamics, where one expects a monotonous 
decay, and for conserved (model B) dynamics, where scaling functions are oscillatory. This is 
achieved through a simple change in the sign of the dimensionful 'mass parameter' //. Although 
we think it unlikely that our non-local representations of ogc(l) should be directly applicable 
to physical models, we consider this qualitative feature encouraging. 

Appendix. On Galilei-transformations 

For comparison with the non-local generators treated in the main text, we recall the computa- 
tion of finite Galilei-transformations, for distributions with a non-vanishing mass ju ^ 0. The 
infinitesimal generator is Y\ii = —td r — fir, from which the finite transformation is formally 
obtained as F(e,t,r) = e~ eYl / 2 F(0,t,r). It is given by the differential equation 

(d E -td r - /ir)F(e,t,r) = , F(0, t, r) = <j>(t, r) (A.l) 

where <fi = <ft(t, r) denotes the given initial distribution. Eq. flA.ip is solved in Fourier space: 

4>(t, k) i — y F(e, t, k) = 4>(t, k + ifj,e) exp 

In direct space, the galilei-transformed distribution becomes 

0(t, r) h-> F{e, t, r) = -L= [ dk 0(t, k + i/ie) e ik(r+te) e^ 12 

= — fdr' <f)(t, r) e ^'e-tite 2 /2 I dk e ifc(r-r'+te) 

v v ' 

2nS(r+te-r') 

= (j)(t,r + te)e' l{r+t£)£ - flt£2/2 (A.3) 

Hence, for fi = 0, the initial distribution is rigidly shifted according to t i — >• t and r i— y r + te. 
This is a consequence of the local nature of the standard Galilei-transformation, which can be 
expressed in terms of a vector field. In particular, the entries in table [1] are recovered. 
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